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Abstract: A simple theory of supersymmetric dark matter (DM) natm'ally linked to 
neutrino flavour physics is studied. The DM sector comprises a spectrum of mixed Ihd- 
rhd sneutrino states where both the sneutrino flavour structure and mass splittings are 
determined by the associated neutrino masses and mixings. Prospects for indirect detec- 
tion from solar capture are good due to a large sneutrino-nucleon cross-section afforded 
by the inelastic splitting (solar capture limits exclude an explanation of DAMA/LIBRA). 
We find parameter regions where all heavier states will have decayed, leaving only one 
flavour mixture of sneutrino as the candidate DM. Such regions have a unique 'smoking 
gun' signature — sneutrino annihilation in the Sun produces a pair of neutrino mass eigen- 
states free from vacuum oscillations, with the potential for detection at neutrino telescopes 
through the observation of a hard spectrum of t'^ and (for a normal neutrino hierarchy). 
Next generation direct detection experiments can explore much of the parameter space 
through both elastic and inelastic scattering. We show in detail that the observed neutrino 
masses and mixings can arise as a consequence of supersymmetry breaking effects in the 
sneutrino DM sector, consistent with all experimental constraints. 
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1. Introduction 

The substantial observational evidence for non-baryonic dark matter (DM) is one of the 
clearest calls for new physics beyond-the-Standard- Model, specifically the existence of an 
exotic stable or very-long-lived massive particle or particles. All this evidence, however, 
is indirect in form, and to-datc docs not fix the nature of the DM particle beyond requir- 
ing it to be colour and electromagnetically neutral and relatively weakly self-interacting. 
Despite this ignorance, one particular candidate has tended, until recently, to dominate 
theoretical speculation — ^the neutralino of the minimal supersymmetric (SUSY) extension 
of the Standard Model, the MSSM. Although the neutralino is the lightest supersymmet- 
ric particle (LSP) over substantial regions of parameter space, and thus stable if i?-parity 
conservation is assumed, the increasing reach of both direct DM searches and collider con- 
straints on the MSSM now pushes MSSM neutralino DM into uncomfortable corners, as 
has been widely recognized. The problems are two-fold: First, direct search experiments 
probe WIMP-nucleon cross sections orders of magnitude smaller than canonical weak scale 
cross sections; and second, given the direct search and collider constraints, the standard 
freeze-out mechanism of dark-matter-genesis generically generates a substantial excess of 
MSSM neutralino DM, limiting the allowed neutralino parameter space to special 'fine- 
tuned' regions of parameter space with coannihilations or resonances. Overall, neutralino 
DM matter works much less naturally than it did in the early 1990's after the excitement 
of the LEP-I results pointing to SUSY unification. 

An independent reason to possibly doubt the standard MSSM neutralino story is that 
it is based on an assumption of simplicity of the dark sector that is probably not warranted 
given our experience of normal matter. In the observable sector there is an extraordinary 
richness of stable or very long-lived massive states, including more than one hundred essen- 
tially stable nuclear isotopes (on the timescale of the Hubble time), a stable charged lepton, 
and three essentially stable neutrino mass eigenstates of various flavour compositions. This 
is despite the fact that all global discrete and continuous symmetries such as individual lep- 
ton number, as well as total lepton and baryon numbers are violated (by at least neutrino 
mixing plus the electroweak anomaly, and likely GUT-suppressed interactions too which 
violate B — L also) . There is no reason why we should expect the DM sector to be any less 
complicated than the observed world, and in a sense the observation of massive neutrinos 
already directly supports this hypothesis — the WIMP dark sector is already composed of 
at least a four-component cocktail of the three massive neutrino species plus one more 
WIMP. 

Motivated by these arguments we, in this paper, investigate the physics of a very slight 
modification of the usual SUSY LSP hypothesis that naturally and simply accommodates 
a much more structured and varied DM sector (and in fact one that has attractive and 
potentially testable connections to neutrino flavour physics). Specifically we study what we 
consider to be one of the simplest alternatives to standard SUSY neutralino DM, though 
maintaining the SUSY framework that naturally possesses the well-known successes of 
weak-scale SUSY theories, such as precision gauge coupling unification and a solution to 
the hierarchy problem, namely, mixed Ihd-rhd sneutrinos. 
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Sneutrinos are potentially interesting alternate SUSY DM candidates since, unlike 
Majorana neutralinos, they can carry flavour quantum numbers while being charge and 
colour neutral, In addition they are the LSP in reasonably large regions of parameter 
space. ^ The traditional reason why sneutrinos have not been considered a good DM 
candidate is that pure Ihd, ie, S'{7(2)2,-active, sneutrinos have too large an annihilation 
cross section, and thus too small a freeze-out relic density to be the observed DM. However 
this problem is easily solved if one posits the existence of weak-scale rhd, ie, SU (2)i-sterile, 
sneutrinos which mix with the Ihd states once electroweak symmetry and SUSY is broken. 
Because of the observed family replication it is most natural to posit three rhd sneutrinos, 
one associated with each lepton flavour, in addition to the three standard Ihd sneutrinos. 

Therefore we are led to a model with six weak-scale sneutrinos, that mix with each 
other and form, as we will argue in detail, a structured DM sector carrying lepton flavour 
quantum numbers identical to that carried by the neutrino mass eigenstates observed in 
neutrino oscillation experiments. In addition, each of these complex scalars splits into 
CP-odd and CP-even components, so the full spectrum of DM states is thus twelve real 
degrees of freedom with a variety of masses. 

Shortening the story of the subsequent sections, what we flnd is a successful model 
comprising an entire sector of mixed sneutrino DM with the novel feature that the sneu- 
trino flavour structure is identical to that of the neutrinos. Further to this the sneutrinos 
have mass splittings between scalar and pseudoscalar components that are proportional 
to the neutrino masses. (See Figure 1.) In addition to these interesting features we find 
regions of parameter space where the splittings are large enough that all states will have 
decayed down, leaving only one flavour of sneutrino as the candidate DM. These regions 
of parameter space potentially have a unique 'smoking gun' signature as pair annihilations 
of these sneutrinos through neutralino exchange produce just one neutrino mass eigenstate 
f3 (for the normal neutrino mass hierarchy), free from vacuum oscillations, that has the 
potential for detection at future neutrino telescopes through the observation of a hard spec- 
trum of z/^ and I'r but not I'e- These neutrinos would have energy equal to the sneutrino 
mass and could arise from annihilations in the Sun. 

Since we require rhd sneutrinos at the weak scale our model also possesses sterile 
rhd neutrinos at this scale or below. Naively one may therefore think that, as the rhd 
neutrinos are not superheavy, it is impossible to generate, eg, by the see-saw mechanism, 
an acceptable spectrum of light neutrinos as observed in neutrino oscillation experiments. 
However as originally show by Arkani-Hamed etal [15, 16] and Borzumati etal [17, 18] 
and amplifled on and extended by later authors (see, eg, [19, 20, 21, 22, 23]) this is not 
the case. In fact there occurs a new and extremely attractive mechanism of neutrino 
mass generation linked to supersymmetry breaking, in a sense generalizing the Giudice- 
Masiero mechanism for generating the /x-term — the Higgsino mass term — in the MSSM. 

^Another attractive feature of considering such sneutrino DM is that it might afford a dynamical ex- 
planation of the baryon-to-dark-matter ratio ^t/^DM — 1/5 as the DM can now possess a lepton num- 
ber asymmetry connected to the baryon-number-asymmetry, in contradistinction to Majorana neutralino 
DM matter which can carry no such asymmetry. We will not pursue these ideas in this work, but see 
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] for studies along these lines. 
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Figure 1: The spectrum of sneutrino masses after lifting of scalar and pseudoscalar mass degener- 
acy by the mass term M^. The red, yellow and blue lines represent, respectively, the electron, muon 
and tau lepton number flavour components, which mirrors the neutrino flavour structure (here we 
have assumed the normal neutrino mass hierarchy). Note that the 12 real degrees of freedom split 
into a heavy sector and a light sector of 6 states each, with further fine structure splittings among 
these two groups, and the splittings between scalar and pseudoscalar components are proportional 
to the neutrino masses. This diagram is not to scale. 

This model of neutrino masses allows for new explanations of the origin of the neutrino 
mass and flavour structure [15, 16, 21, 23] and elegantly accommodates such nice features 
as weak-scale resonant leptogenesis [22]. Since our main focus in this work is the novel 
sneutrino DM phenomenology we only very quickly summarize the physics of neutrino 
mass generation in as far as it impacts the sneutrino sector, and we strongly encourage the 
reader to refer to these other papers for more details. 

Before we start, one other feature that deserves discussion is the fact that this class 
of models was the original and motivating example of inelastic DM [19, 20], providing 
a possible explanation of the reported DAMA/LIBRA [24] signals consistent with the 
exclusions reported by other direct DM detection experiments. Since these original models 
involving sneutrinos, the idea of inelastic DM has been implemented in many more set ups, 
and has been much studied (see, eg, [25] and references therein). This inelastic scattering 
also changes the phenomenology of DM solar capture and we include the exclusion limits 
from Refs. [26, 27] on the nucleon scattering cross-section. These can be more constraining 
than the current generation of direct detection experiments over large regions of parameter 
space. 

There are many previous studies of sneutrino DM, see eg, [28, 29, 30, 31, 32, 33, 34, 
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35, 36, 37, 38, 39, 40, 41, 42, 43, 44], though most do not consider connecting the sneutrino 
and neutrino flavour structure. (For treatments of sneutrino DM taking hmits from the 
generation of neutrino masses we refer the reader to [33, 34, 44].) 

The previous work most closely related to that considered here is Ref.[44] where the 
flavour structure of the sneutrinos is also considered. However our model has the novel 
feature that the neutrino and sneutrino flavour structures are exactly the same, allowing a 
predictive 'smoking gun' signature for sneutrino annihilations into neutrinos, as discussed 
in Section 4.1. Another important difference is that in our model there is an interplay 
between two terms contributing to the neutrino mass generation (see the off-diagonal entries 
of eqs.(2.4)) that allow us to fit the observed neutrino masses consistent with successful 
sneutrino DM over a significant region of parameter space. Finally, our analysis includes 
the important limits arising from solar capture, which excludes the possibility of an inelastic 
DM explanation of DAMA/LIBRA using mixed sneutrinos. 

In Section 2 we first discuss the specific model, outlining the generation of neutrino 
masses and the origin of the identical neutrino and sneutrino flavour structure, as well 
as the resulting mass spectrum of the sneutrino states. We then go on in Section 3 to 
present the decay lifetimes for the different flavoured sneutrino states and consider the 
DM phenomenology of the sneutrinos, showing the regions of parameter space where they 
are good DM candidates. In Section 4 the potential unique signatures of this model are 
considered, while our conclusions are contained in Section 5. Two short appendices contain 
technical details. 

2. The Model 

We start by summarizing the relevant fleld content and effective Lagrangian of the model 
we employ — that of [21], to which we refer the reader for additional details. 

2.1 Field content and interactions 

To the fleld content of the MSSM we append three SU{3) x SU{2) l x C/(l)y-sterile neutrinos 
and their sneutrino superpartners, hi, which combine into the chiral superflelds Ni. Here 
i = 1,2,3 labels the generation. The terms arising from the superpotential are 



The SUSY-preserving superpotential term for the rhd neutrino mass, M]\fNiNi, arises 
from a higher-dimension Kahler term involving a SUSY-breaking spurion F-component, 
in analogy with the Giudice-Masiero mechanism for the //-term, thus giving rhd neutrino 
masses at the weak scale Mat ~ rnj/M ~ TeV. (Here m/ ~ 10^*^ — 10"*^^ GeV is the 
intermediate scale at which SUSY is broken in the hidden sector, while M is the reduced 
Planck scale.) The Yukawa couphng Ajj between the Ihd lepton doublet superflelds, Li, and 
the rhd neutrino supermultiplets, Ni, also arises from a Kahler term involving a spurion 
F-component and is suppressed in magnitude by a factor of \Xij\ ~ rrii/M ~ 10""^ — 10~^. 
These suppressions can be justifled with an R-symmetry [21]. 




(2.1) 
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In addition to the usual soft-SUSY breaking terms of the MSSM there are also SUSY- 
breaking terms of the form 



A£ = m| |nj|^ + ALihihu + -Ay M'^hihj + h.c. . (2.2) 



Namely, TeV-scale, but flavour diagonal, soft mass and trilinear scalar A-terms, and a 
small but significant rhd sneutrino lepton-number violating B-term with coefficient Bfj ~ 
XijM^ ~ m]{mi/MY ~ (TeV^/M)!/^ ~ (100 MeV)^ and flavour structure identical to 
that of the neutrino Yukawa coupling, which we assume to be real. Again this structure of 
suppressions and lepton flavour breaking can be justified by the R-symmetry and flavour 
structure properties of the SUSY- breaking spurious. More complicated flavour breaking 
patterns, for instance those with non-diagonal ^-terms, are also possible, but eq.(2.1) and 
(2.2) is the simplest structure that leads to a successful spectrum of masses and mixings 
in the neutrino sector (as we quickly recall in Section 2.3). It also leads to a simple and 
direct connection between the neutrino and sneutrino flavour structure and masses. 

2.2 Sneutrino Masses 

In total we have six complex sneutrino fields, however these are split into scalar and pseu- 
doscalar components by the sneutrino i?-term in eq.(2.2), leaving twelve mass eigenstates. 
After electroweak symmetry breaking the sneutrino mass matrix has the form: 



( Mlly, At;sin/?l3 \ijMNVS\n(3\ 

Avsin/3l3 M|l3 \ijMNVsm.[3 \jMl 

\ XijM^v sin 



(2.3) 



where the basis is (i>*, n, u,fi*). Due to the flavour diagonality of the unsuppressed terms 
after electroweak symmetry breaking, this 12 x 12 sneutrino mass matrix is made up of 
a total of sixteen 3x3 blocks which are either proportional to the identity matrix (from 
the A-term, Majorana mass F-term and D-terms), or Ajj (from the B-term and the rhd 
sneutrino F-term). Therefore the sneutrino flavour structure is completely determined by 
the matrix that diagonalises Ajj, which we will call Ux. 

After diagonalisation of the flavour structure the sneutrino mass matrix is in the form 
of sixteen diagonal 3x3 blocks, which by further rotations can be manipulated (detailed 
in Appendix A) into the simple block-diagonal form of six 2x2 matrices, three for the 
scalars and pseudoscalars each: 

=( Av sin f3 + X^a,M],\ 

yAvsmp + Xaa^Ml M^ + X^M^ J ^ ' ' 

=( AvsinP-X^a,Ml\ 

yAvsmfi-Xa^ayMl M\-Xa,Ml j ^ ' ' 

where v'^ and v~ denote the scalar and pseudoscalar components respectively. The sub- 
script a = {1,2,3} denotes the mixed flavour sneutrino mass eigenstate (note the change 
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from Latin to Greek indices on going from weak to mass eigenstates), Aq, are the eigenval- 
ues of the Xij matrix, M| = m| + cos 2/3 and = m| + M^, where rrij^ and 
are the usual soft mass terms, and, v = 174 GeV, is the Higgs expectation value. We have 
also defined 

vMNsmf] 

a. = (2.6) 

which gives a measure of the relative size of the B-term and F-term contributions to the 
splitting of the CP-eigenstate masses. 

These 2x2 matrices can then be diagonalised to complete the exact sneutrino rotation 
matrix U . Further, in this form one can see there are twelve mass eigenstates: two z/+ and 
two u~ for each flavour, and both the heavy and light and u~ states for each flavour 
have their masses split due to Aij, becoming degenerate in the limit A^ — > 0. Expanding to 
first order in the small parameter A^ the splitting between the light v'^ and vj^ components 
for each flavour is given by: 

5a,L = M-+ - M-_ = A„— ^(cos^ e - a^sm29) ~ 100 keV (2.7) 

' oc,L a,L Myj^ 

where My^ and 6 are the lightest mass eigenstate and the rotation angle which diagonalises 
the 2x2 matrices in eq.(2.4) in the limit Aq, — > 0. Similarly, the splitting between the and 

vjj components is given by 5a^H = Af,>+ — = Xa-^^{s>va^ + a^sin2^) ~ 10 keV, 

where Mj>^ is the heavy mass eigenstate in the limit Aq, — > 0. Therefore we are left with a 
sneutrino mass spectrum as shown schematically in Figure 1. 



2.3 Neutrino Masses 

As the Majorana mass for the rhd neutrinos (M^r ~ 1 TeV) is very small compared to the 
standard seesaw set-up (Mat ~ 10^*^ — 10^^ GeV) one would naturally expect the seesaw 
contribution to the light neutrino masses to be large, however here this contribution is 
suppressed by the square of the small neutrino Yukawa coupling (A^ ~ 10^^^ — 10^^^) 
and instead the dominant Majorana masses for the neutrinos arise radiatively through 
loops involving sneutrinos and neutralinos. This is due to the combination of the lifting 
of the degeneracy between and v~ sneutrino states and the Majorana nature of the 
neutralinos and can be seen schematically through loop diagrams involving perturbative 
mass-insertions, as in [21]. However, as it is possible to solve for the sneutrino masses and 
mixing matrices analytically, the exact one-loop contribution to the neutrino masses can 
be calculated from the Feynman diagram in Figure 2. 
Correspondingly, the neutrino masses are given by; 

1 /M \^ ^ 

^"^^ " 2 (477^) Y.Y.'^xMxi^l^ew - N^2Cosew?ul^L{r^,rn^^)U^,,+Q (2.8) 

^ ^ a=l x=l 
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) = + 1 + 1^ ( ) (2-9) 



/ \ 

I'i * ' » < 1 < 
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Figure 2: Feynman diagram corresponding to the radiative neutrino mass generation. 

where a runs over the twelve sneutrino mass eigenstates, N^y are the neutralino mixing 
matrices, the neutralino masses, = Ma/m^^ and:^ 

The parameter dependence and flavour structure of the neutrino masses is more ap- 
parent by expanding the previous exact one-loop formula to first order in the \ij matrix. 
This gives; 

Xij ( — - ) (iV^i sin 9w - Nr,2 cos 9wf 

V 47ri' y ^ m^, (2.10) 

• (sin^ 0(cos^ — ttu sin 29)Li{xx,yx) + cos^ 6'(sin^ 6 + a,^ sin 29)L2{xrc, Ux)) 

where Li and L2, defined in Appendix B, are the loop contributions from the sneutrino 
soft mass term and the F-term contributions. The parameter a^, defined in eq.(2.6), gives 
the freedom to generate neutrino masses that aren't too large by canceling these two 
contributions off each other. 

From eq.(2.10) one can see that the neutrino masses are proportional to the eigenvalues 
of Xij. Combining this with eq.(2.7) it is clear thar 6a oc Mjy^. In order to obtain the correct 
neutrino masses it is necessary that \a/\p = M^^/M^^ and therefore 6a/ 6/3 = M^^/M^^. 
Once these ratios and all other parameters are set, then the overall magnitude of the 
neutrino masses depends on the parameter a^. We allow to take values — 1 < < 1, 
and find that for reasonable values of this parameter, neutrino masses below cosmological 
bounds and with M^^ > Y^Amf^~+~Am|^ are obtained. This will be discussed in more 
detail in Section 3.3. 

It is clear from eq.(2.10), and is shown to all orders in Ay in Appendix B, that the 
flavour structure of the neutrino mass matrix takes a simple form. The PMNS neutrino 
mixing matrix is given simply by Upmns = ^Ai where U\ was defined previously as the 
matrix which diagonalizes Xij , and also completely determines the flavour structure of the 
sneutrinos. Therefore the flavour structure of the neutrinos is automatically identical to 
that of the sneutrinos in this model. This enables us to make definite statements about 
possible 'smoking gun' signatures. 

Hence, in order to reproduce the measured neutrino mixing and mass parameters the 
form of Xij is completely known, up to an overall normalisation factor which is absorbed into 



^As X^Q^i Ua,iU^ j^a = Sij the last two terms vanish when summing over the sneutrinos as they are 
independent of the sneutrino masses and a. 
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M^. Since 6a oc M,^^ , the largest sneutrino splitting 6a is in correspondence with the largest 
neutrino mass. For the normal neutrino mass hierarchy, we know the heaviest neutrino's 
flavour structure is almost exactly half fj, and r, therefore we know the four sneutrino states 
i/3 (two heavy scalar and pseudoscalar, and two light scalar and pseudoscalar states) must 
also have this flavour structure. This is shown schematically in Figure 1. 

3. Dark Matter 

In the previous section we have derived the sneutrino mass spectrum and shown that, 
as well as giving mass to the neutrinos, they also pass on their flavour structure to the 
neutrinos. In this section we would like to apply current experimental constraints to flnd 
the regions of parameter space where the flavoured sneutrinos are a good DM candidate.^ 
As explained above, the sneutrino spectrum is split into two groups of six particles with 
the same flavour structure as the neutrinos. In this section we will only be concerned with 
the possible transitions of the six lightest particles as the heavier ones very rapidly decay. 
We will find that there are three general possibilities for the current DM composition: 1) 
In a Hubble time all of the heavier sneutrinos could have decayed to the lightest state, 
^s'l' which makes up all of the measured DM relic density; 2) The three lightest species 
survive to the present day, each making up a third of the measured relic density; 3) An 
intermediate situation where one or more of the heavier sneutrino states have decayed (or 
are currently decaying) down to the lighter states, but a non-trivial cocktail of different 
flavoured sneutrinos remains. For pedagogical clarity, we shall assume in the following 
that the neutrinos follow the normal mass hierarchy, however we will briefly highlight the 
differences if the inverted mass hierarchy is assumed at the end of this section. 

3.1 Stable and meta-stable sneutrino spectrum and decays 

The lifetime of the lighter six states is a sensitive function of the splittings 61 and Ai, 
(deflned in Figure 3), which as shown in eq.(2.7) depend on the ratio a^. As we will 
show later, there are regions of the parameter space of Qi, where the lifetime of the decays 
of D^j^ and to z/g"^ are longer than the age of the Universe. In this case, the DM is 
comprised of equal abundances of all three states t']"^, 1^2 L ^^'^ ^3L^ ^^'^ therefore overall 
flavour neutral. However there are also portions of parameter space where the states will 
have decayed to only t'^j^, leading to one flavour of sneutrino as the DM candidate with 
potential signatures linked to the mixing properties of neutrinos. 

Initially we consider the decay of the lightest scalar particle, i^ii, as shown in Figure 3. 
The particle can decay by emission through intermediate Z exchange to the pseudoscalar 
particle of the same flavour, i^f^, or by z^z^ emission through an intermediate neutralino 
exchange to any of the lighter pseudoscalar particles. Typically the dominant decay mode 
is through the Z, 



^For all of our studies we take the current total DM density to be given by the WMAP5+BAO+SN 
value nh'^ = 0.1143 ± 0.0034[45]. 




(3.1) 
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Figure 3: The six lightest sneutrino mass eigenstates (sphttings not to scale). can decay 

to iy^j^ by vv emission through intermediate Z exchange, as shown by the blue line, or to any of 
the lighter states by w emission through a neutralino (red lines) . and 1^2 l can decay to v^j^ 
through an intermediate neutralino exchange. For some regions of parameter space, these decays 
will be longer than the age of the Universe. 

and is around an order of magnitude faster than neutralino-mediated decay. The associated 
lifetime is, for all reasonable parameters, short enough that no i^^^ states survive to 
the current day. Although these decays occur after Big Bang Nucleosynthesis there are no 
changes in the predicted light element abundances as the mass splittings are small and there 
is no hadronic energy injection. There are similarly no constraints from spectral distortions 
in the CMB [46] from the energy injection due to Z-mediated decays to two photons, as 
the loop-suppressed branching fraction into photons is too small for these features to be 
observable. 

We will now consider the decays of the pseudoscalar particles, and in particular, the 
decay of z/j~^ to P^^, as indicated by the dotted red arrow in Figure 3, since it will be the 
fastest of the possible decays. If this lifetime is longer than the age of the universe, the 
lifetime for i>^^ to decay to D^j^ will also be longer. These decays can only occur by 
emission through an intermediate neutralino. The decay width for typical parameters is: 

~ (10- yrs)- (^) f f— V r^^V f ?2^V (3.3) 

where x = (Ex ^i^xi sm9w - N,2 cos ^h/)')'+(Ex ^^^l^xi smOw - N^2 cos 9w\^) 
depends on the neutralino mixing matrix, N^y, and neutralino masses, m^^. 

The decay width eq.(3.3) is in the range where, depending upon parameters, the pseu- 
doscalar states can either have essentially all decayed by the present epoch, or can 
be substantially still present, or can be currently decaying. For the purposes of matching 
to the observed DM density, however, this ambiguity is immaterial as each heavier state 
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decays to one only marginally lighter state. The same is true for the decays of the ^ 
states, so in computing the total freeze-out relic density it is a good approximation to sum 
the contribution of all six light states ignoring the later effect of decays. 

3.2 DM relic density and experimental constraints 



Although various thermal and non-thermal mechanisms for generating the observed 
DM relic density are possible, we assume that the relic density is generated by the standard 
thermal freeze-out process. In particular we do not here consider the possibility of using the 
calculable and IR-dominated thermal freeze-in process recently advocated by Hall etal [47] , 
as this would apply for supersymmetric theories with Dirac neutrino masses, rather than 
the model outlined in Section 2. 

The three main processes which contribute to the freeze-out relic density calculation 
are shown in Figure 4. We performed the calculation using micrOMEGAs 2.2 [48] with 
model files created using LanHEP 3.0.4 [49] and the results are shown as the solid green 
curves in Figures 5, 6 and 7. Coannihilation between the six lightest sneutrinos, which is 
obviously important since they are almost degenerate in mass, is automatically included 
in micrOMEGAs. The dashed black and purple lines show the limits from the current 
generation of direct detection experiments on the DM-nucleon scattering cross section, or 
equivalently from using eqs.(3.4) and (3.5), on the mixing angle sin0. An interesting feature 
of this model is the possibility of two signals at future direct detection experiments; the DM 
can up-scatter inelastically via Z exchange, for example from to and therefore is 
sensitive to the small mass splitting 5a,L^ or can elastically scatter by Higgs exchange. This 
also means that we can constrain the parameter space in two regions due to the different 
kinematics of each collision. The elastic cross section is larger at smaller masses Mp^ , while 
the inelastic cross section limit is stronger for smaller 5a,L, as we will explicitly demonstrate 
later. The cross section for coherent scattering off a nucleus at zero momentum transfer 
by Z exchange is 



Here /i is the reduced mass for the sneutrino-nucleus, Z the proton number and the 
number of protons and neutrons in the nucleus. The cross section for elastic scattering off 




Figure 4: The dominant processes which set the sneutrino relic density. 




((yl-Z)-(l-4sin2 ew)zf 



(3.4) 
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Figure 5: Exclusion curves and thermal relic density constraints in the sin 6* - Af^^ plane for 
Ml = /i = 250 GcV, Ma = 500 GeV and tan/3 = 10, and with Higgs mass 121 GeV. The green 
solid line indicates the value of sin 9 where the observed DM relic density is generated thermally. 
The dashed lines give the exclusion curves from direct detection experiments and the shaded regions 
show the values of sin which fit the DAMA/LIBRA experiment at 90% and 99% confidence levels. 
The dot-dashed curve gives the exclusion curve from the measurement of solar neutrinos from Super- 
Kamiokande. The green solid line below the dashed and dash-dotted lines indicate parameters where 
the sneutrinos are a good thermal DM candidate. The two features below 60 GeV on this thermal 
relic density curve correspond to the Higgs- and Z-funnel regions. The Z-funnel is a co-annihilation 
effect and so would be absent from indirect detection signals. Note that, with the exception of the 
Higgs funnel region, the elastic direct detection, and solar capture indirect detection limits are close 
to the region of parameter space required for a successful thermally generated relic density. 



a nucleoli by Higgs exchange in the decoupling limit is given by: 

aivnf VfAsin/3sin2g-^cos2/?sin^gV (3-5) 

Sirml ynin + rrij^ J \ v J 

where m„ is the nucleon mass and ghnn = (YI^''^'^ It^ + 3/ YI'q'''^ Itg^ ■ We use the 
values of /t^ and Jtq from Ref.[33] which give ghnn = 1-43 x 10^'^. 

The red dot-dashed lines show the indirect detection limits on sin 9 from the DM 
capture, and subsequent annihilation into neutrinos, by the Sun. At the masses we consider, 
the limits are set from the observation of solar neutrinos by Super-Kamiokande [50]. Note 
that the solar capture and annihilation limits are the most constraining limits over a sizeable 
region of parameter space. 

In the limit Aq, — > 0, the sneutrino masses and mixings are determined by Ml, Mr, A 
and tan /?, however we find it more intuitive to trade in Ml, Mli and A for M^, Mi,^ and 
9. Furthermore, once and one of the 5a,L are specified, and Ml, M2 and (i are chosen 
so that the neutralino masses m^^ and mixing angles N^y can be computed, the absolute 
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neutrino masses can be calculated from eq.(2.10) and the other two 5a, l can be found from 
the ratio 6a,L/5/3,L = ^val^vp- -^t the top of each figure we have displayed the value of 
and 6 we have fixed. Note that 5 always refers to the smallest (5a, l from the sneutrinos 
which have not decayed. 
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Figure 6: Left panel: Exclusion curves and thermal relic density constraints for Mi = /i = 
150 GeV, Ah = 300 GeV and tan/3 10. As well as the Higgs-funnel region at - 60 GeV, the 
Z- funnel region at ^ 45 GeV is now also viable. Again, note that, with the exception of the Higgs 
funnel region, the elastic direct detection, and solar capture indirect detection limits are close to the 
region of parameter space required for a successful thermally generated relic density. Right panel: 
The branching fractions for the same parameters as in the left panel. For each value of A'/p^ , we 
have chosen sin^? to give the correct thermal relic density. At low masses the branching fractions 
to neutrinos dominates, and at masses just below the W mass is substantial. At masses above the 
W mass, that to W~^W~ dominates. 

For all further calculations we fix the MSSM parameters at the weak scale. We set 
tan/3 = 10, the pseudoscalar Higgs mass Ma = 500 GeV, the rhd slepton soft mass = 
250 GeV, the top soft coupling At = —1 TeV and all other soft parameters to 1 TeV giving 
a Higgs mass of 121 GeV. Ml, M2 and n are varied in different plots but we maintain the 
relation M2 = 2M1. 

The dominant channel which contributes to the relic density is through the Higgs s- 
channel because of the large A-terms. An especially noticeable feature in the relic density 
curve in Figures 5, 6 and 7 is the Higgs funnel at half the Higgs mass ~ 60 GeV. A similar 
effect can be seen at ~ 45 GeV, the Z funnel, but is much smaller because it is suppressed by 
sin^ 6 rather than sin^ 2^, cf. eqs.(3.4) and (3.5). The aimihilation through the neutralino 
t-channel is important at lower mass (< 45 GeV) and just below the W mass, where, from 
the right panel of Figure 6, we see that the branching fraction into neutrinos can be large. 

As increases with 6 and Mp^ held constant, the A-term increases, so smaller values 
of sin 9 are required to ensure the sneutrinos do not over-annihilate producing a relic density 



- 12 - 



1.00 



0.50 



0.20 



0.10 



0.05 



0.02 



6 = 50 keV, mi = 400 GeV 



I.OOr 



6= 150 keV, ni^ = 400GeV 




0.1109 < rih- < 0.1177 
Direct Detection (Inelastic) 
Direct Detection (Elastic) 
Solar Capture 
DAMA/LIBRA Region 



100 150 
Mf, [GeV] 



100 150 
Mj,^ [GeVl 



200 



Figure 7: The effect of changing 5 while keeping everything else fixed. The rehc density and 
elastic direct detection curves remain the same, while there is a small change at low masses in 
the solar capture limits. The major change is in the inelastic direct detection limits which change 
considerably due to the different kinematics of the collision with the target nucleus. 



below the measured value. From eq.(3.5), we can also see that this increase in A is the 
reason why the elastic direct detection limits become stronger as rrij^ is increased. 

For < 350 GeV, the only region of parameter space not excluded from elastic direct 
detection and solar limits is below the W mass; as well as the Higgs- funnel region, the left 
panel of Figure 6 demonstrates that the Z-funnel region at ~ 45 GeV also opens up. For 

> 500 GeV, the elastic direct detection limits and relic density curve move to lower 
values of sin together so that no new lower values of m^;^ are allowed. 

Changing S has virtually no effect on the relic density calculation and the elastic 
direct detection limits, while the change in the solar limits is very minor. However as 
Figure 7 shows, the major difference is in the inelastic direct detection limits which change 
considerably. 

The somewhat erratic appearance of the solar limits in Figures 5,6,7 and 8 is because 
we have only plotted the strongest constraint arising from the different possible annihila- 
tion channels; i^z/, rr, bb, W'^W^ , ZZ and hh. Typically we found that for M;>^ < Mw, 
the limits from annihilation into I'u dominate except for when Mj>^ ~ Mh/2, where the rr 
limits dominate. This is because the neutrinos are produced through t-channel neutralino 
exchange which are not enhanced in the Higgs funnel, as other massive particles are. For 
Mi)j^ > M\y, the W~^W~ limits typically dominate. The branching fractions were calcu- 
lated using micrOMEGAs and we used the limits on the annihilation branching fraction 
derived using the methods in [26]. We have also only shown the strongest limits from the 
direct detection experiments. We assume the standard halo model and refer the reader to 
Ref. [51] for details on how these limits were calculated and the experimental details. 
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Figure 8: The grey shaded regions show the values of the scattering cross section a„ which fit 
the DAMA/LIBRA experiment at 90% and 99% confidence levels. Also shown are the best direct 
detection limits (purple dashed) and solar capture limits (red dot dashed). Clearly the solar capture 
limits completely exclude the DAMA/LIBRA region in this model. 

Since our model naturally includes a splitting of order 100 keV, it is natural to ask 
whether this is consistent with the inelastic DM [19] explanation of the positive annual 
modulation signal measured by the DAMA/LIBRA [24] experiment and the null results 
from the other direct detection experiments (notably CDMS II [52] and XENONIO [53]). 
Although DAMA/LIBRA may still be consistent with the other direct detection exper- 
iments at 5 ~ 125 keV [54, 51, 25, 55, 56], we find that when the solar capture limits 
from Refs. [26, 27] are included as in Figure 8, the DAMA/LIBRA region is conclusively 
excluded in this model. 

3.3 The allowed Cj/ parameter space 

For a given set of the parameters {m^^, /i, Ml, M2, tan /3} one can see from Figures 5, 6 
and 7 that the allowed values of sin 6 and M^j^ are tightly constrained by limits from direct 
and indirect detection. Although non-thermal mechanisms for generating a relic density 
are possible, we also impose that the correct relic density is generated thermally and this 
reduces the allowed parameter space further. Therefore once is chosen, sin 6 is set by 
this constraint. 

Finally it remains to constrain the parameter which contributes to the inelastic 
splittings and the neutrino masses through eqs. 2.7 and 2.10. For a given splitting ^3^^ 
the allowed parameter space for a^, is determined by the requirement that the neutrino 
masses satisfy both the cosmological bound, X^m^ < 0.61 eV [57], and the mass bound 
M^3 > Y^Amf2 + Am^g. 

In Figure 9 we show the allowed values of for a given parameter set with varying 
WIMP mass M/>^ . For a typical set of soft mass and electroweak parameters it is in general 
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M-„^ [GeV] 

Figure 9: Allowed regions of parameter space for the a,^ parameter for /it — Ml — 250GeV, 
M2 = SOOGeV, mp^ = 350 GeV, sin 9 = 0.2, tan /3 = 10 and ^ = 100 keV. 

possible to find a value of a^j for which the neutrino masses satisfy the bounds, although 
the allowed region may become small. 

We have shown that within this model it is possible to generate all of the observed 
neutrino parameters and satisfy cosmological mass-sum bounds. Concurrently the mixed 
sneutrinos have the same flavour structure as the neutrinos and are good DM candidates 
capable of providing the observed relic density and avoiding current direct and indirect 
detection limits. 

It is now interesting to consider properties of the sneutrinos that could lead to a positive 
identification at future experiments. A potentially unambiguous annihilation signal could 
arise if the DM sneutrinos have all decayed down to the lowest state. 

By choosing the same parameter set as that used to generate Figure 5 and varying 
values of svaO and WIMP mass, we have calculated the fraction of the allowed parameter 
space of the parameter ay for which the sneutrinos will all have decayed down to the lightest 
state. The results for two different values of 5 are shown in Figure 10. From this it seems 
that, for this model of DM, the fraction of parameter space for which the DM sneutrinos 
are made up of one state lies in the range 5 — 30% for WIMP masses above Mw and the 
range — 45% for masses below Mi^. 

3.4 Neutrino mass hierarchy 

This model does not have a preference for the neutrino mass hierarchy, since it is equally 
capable of generating either. In the discussion above we have assumed the normal mass 
hierarchy. Here we highlight how the phenomenology changes if the inverted hierarchy is 
assumed. Since Sa^L oc My^, the 6a,L follow the neutrino mass hierarchy so that in the 
normal hierarchy, d-^^i is the largest splitting, and is much bigger than 52,l and which 
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Figure 10: Contours of constant fraction of a^, parameter space for which the sneutrinos have 
decayed down to the lightest state (i.e. constant fraction of allowed values of Ui, for which 
< lO^'^ yrs). For comparison the regions where the observed relic density is produced 
are superimposed. The SUSY-breaking parameters are /x = Ml = 250 GeV, M2 = 500 GeV, 
tan (3 = 10 and the solid and dashed lines are for = 350 and 500 GeV. 



are expected to always be close together because (^^l ~ '^iL °^ ^"^ii) which is known to 
be small from solar neutrino measurements. Conversely in the inverted hierarchy, ^3 ^ will 
be the smallest splitting, 62,1 will be the largest, but 5i^l will be similar because Am^i 
is small. Therefore even if u^j^ has decayed, we expect to have a significant fraction of 
both D^j^ and and the overall flavour structure of the DM will depend precisely on the 
relative abundances of 1^7 r and i>-7r- 



4. Detection and Identification 
4.1 Indirect Detection 

Mixed sneutrinos present an interesting scenario for indirect detection through neutrino 
annihilation products for three reasons: 

• The inelastic splitting allows for a comparatively large WIMP-nucleon scattering 
cross section, due to the kinematic suppression of scattering at direct detection 
experiments. However, the large escape velocity of the Sun implies a typical WIMP 
kinetic energy in the Sun greater than this inelastic splitting, thus allowing a relatively 
large solar capture rate. 

• Sneutrino annihilation can proceed via t-channcl neutralino exchange to two neutrinos 
which propagate relatively unhindered (with some attenuation in the Sun) to the 
Earth. These neutrinos could be detected in future neutrino telescopes as a hard 
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spectrum centered at the sneutrino mass, although this hard spectrum may he below 
the flux of neutrinos from annihilations to other SM particles. This is illustrated in 
Figure 11. 

• If the DM were made up of all three lowest lying states then annihilations to neutrinos 
would result in equal numbers of v^, and v^-, with the directional flux of hard 
equalling roughly one third of the total flux of hard neutrinos. If the sneutrinos have 
all decayed down to the lightest state then the hard spectrum of neutrinos would be 
comprised solely of the 1^3 mass eigenstate, and the flux would equal roughly one 
half of the total flux for the normal hierarchy neutrino structure. In the vacuum none 
of the neutrinos would oscillate as they are produced in the mass eigenstates i'i,2,3- 
Moreover, if annihilations are purely to ^3 states in the sun, MSW effects would be 
negligible due to the small component of 1^3. 

As the WIMP-nucleon scattering cross section is typically large, sneutrino capture and 
annihilation in astrophysical bodies provides the most promising scenario for indirect de- 
tection. However it should be noted that neutrino signals from annihilation in the galactic 
halo are not subject to the large uncertainties relating to the DM structure in the cen- 
ter of the galaxy [58] , and a dedicated analysis by neutrino telescope collaborations could 
place interesting limits on thermal relic WIMPs from annihilations resulting in flnal state 
neutrinos. 

Due to the suppressed capture rate in the Earth from the inelastic splitting, we will 
only discuss potential signals from sneutrino capture in the Sun. Although a thorough 
treatment of limits from upcoming neutrino telescopes is best left until after the collection 
of data, it is interesting to speculate about the detection potential. Compared to existing 
Super-Kamiokande limits [50] projections of the performance of the completed IceCube -|- 
DeepCore detector indicate an order of magnitude increase in sensitivity to i^^ from the 
sun for WIMP masses ^ 100 GeV [59]. An order of magnitude increase in sensitivity to (T„, 
corresponds to an increase in the limits on sin^ by a factor ^ 1.8. As the regions where 
the observed relic density is generated are already close to current indirect detection limits 
then inspection of Figures 5, 6 and 7 would suggest that the completed IceCube detector 
will be sensitive to large portions of parameter space of this model for < 500 GeV. 

In Figure 11 we plot the neutrino spectrum at the Earth for various sneutrino masses 
where current direct and indirect detection bounds are just satisfled. Inspecting Figure 
6 one can see that when My^^ > Mw (respectively Mz) then W (and Z) bosons are the 
dominant annihilation products, while for Mi>^ ~ then hb dominate due to the Higgs 

resonance. Neutrinos from both these annihilations swamp the hard spectrum of neutrinos 
from neutralino exchange. However when the sneutrino masses lie out of these ranges the 
sharp peak from annihilation to neutrinos can be seen where the annihilation to neutrinos 
is at its greatest. Observation of such a feature would constitute a compelling case for 
sneutrino DM as direct annihilation of neutralinos to neutrinos is forbidden due to their 
Majorana nature. As described in Ref. [38], the peak from annihilations to neutrinos would 
lead to a linear component in the muon energy spectrum at IceCube and, if large enough, 
could be distinguished from annihilation to other standard model particles which give a 
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Figure 11: The spectrum of neutrinos from the Sun per sneutrino annihilation for two sets of 
parameters lying just below current indirect detection bounds. The SUSY-breaking parameters 
are as in Figure 6 and is varied in order to stay just below current indirect detection bounds 
for higher sneutrino masses. The neutrino spectrum is evaluated using the branching fractions 
calculated using micrOMEGAs 2.2 [48] along with the oscillation and propagation code publicly 
available [60]. For the two Ihd panels the hard feature corresponds to two-body annihilation and 
gives neutrinos at energy Mq^ . The branching fraction to neutrinos is large away from the Higgs 
funnel and below the W mass, as shown in the rhd panel of Figure 6. Note that these are raw 
spectra. In the rhd panels we illustrate cases where the direct 2-neutrino final state branching 
fraction is small, the spectrum of neutrinos being dominated by secondary decays of annihilation 
products. 



non-linear muon energy spectrum. Further, if the flavour content of this component or 
alternatively just that of the highest energy neutrinos could be determined and was found 
to be non-democratic, then this would be an indication of neutrino-flavoured sneutrino 



- 18 - 



DM. 

We emphasize that if the sneutrinos have ah decayed down to the Ughtest state then 
the prediction of this model for the ratios of the flavour content of the hard neutrino 
spectrum from annihilations through neutralinos is ~ (0, 1/2, 1/2), in the basis {fe, ffi, fr), 
exactly matching the flavour content of the heaviest neutrino in the normal hierarchy. If 
the inverted hierarchy is assumed and the sneutrinos have decayed down to the two lightest 
states, which are very long lived due to their small mass splitting, the flavour content would 
be expected to take the form ~ (1, 1/2, 1/2). 

4.2 Direct detection 

Referring the reader back to Figures 5, 6 and 7 one can see that current direct detection 
bounds are close to large portions of parameter space for this model, whether for elastic 
scattering via Higgs exchange or inelastic scattering when the splitting 6 is smaller than 
~ 100 keV. Therefore the prospects for the direct detection of sneutrino DM are good as 
future direct detection experiments such as EURECA [61], XENONIOO [62] and LUX [63] 
are expected to have significant increases in detection sensitivity. 

Regardless of the size of the inelastic splitting the elastic scattering via Higgs exchange 
is unavoidable and future direct detection experiments should be capable of placing strong 
limits or possibly observing a positive detection signal with this channel. Also, if the in- 
elastic splitting isn't too large (5 ^ 150 keV) then there is also good potential for observing 
inelastic scattering, which could be discriminated against elastic scattering by the shape 
of the event spectrum, thus giving a unique window onto the size of the inelastic splitting 
among sneutrino states. 

Therefore upcoming direct and indirect detection experiments will provide a compli- 
mentary test of the validity of this model of DM. 

4.3 Collider Signatures 

Due to the small mass splittings of the six lightest sneutrino states, to a good approximation 
all six states would be produced in equal multiplicity at the LHC, making a determination 
of flavour structures at the LHC very difficult. However it may be possible to identify sneu- 
trinos as the LSP. The LHC signatures of mixed sneutrinos have been studied previously 
[34] where it was found that mixed sneutrinos can be distinguished from the MSSM over 
large regions of parameter space. If the NLSP is a right-handed slepton then large lepton 
multiplicities may result from the decays of more massive sparticles. In particular decays 
of the lightest neutralino lead to opposite-sign, same-flavour dilepton signatures. This can 
be distinguished from similar MSSM scenarios by observing the shape of the dilepton in- 
variant mass distribution [34] . For large sin 9 decay chains starting with squarks can lead 
to jet-lepton signatures which are present in the MSSM, but could be distinguished by 
observing sneutrino-charged slepton mass splittings which arise only through electroweak 
symmetry breaking in the MSSM [34] . 

Finally, if the sneutrinos were to lie in the Z-funnel region, M^;^ ~ 45 GeV, as in the left 
panel of Figure 6, Higgs searches at the LHC and ILC could be dramatically altered because 
the large A-term implies that the dominant decay of the Higgs will be to two sneutrinos. 
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For the parameters of Figure 6 and sinO = 0.2, the Higgs decay width Th will increase by an 
order of magnitude from F/i ~ 5 MeV to Th ~ 50 MeV. Although not directly measureable 
at the ILC, this could be determined indirectly from = T{H WW*)/Br{H WW*) 
where Br{H — > VFM^*) is directly measured and T{H — > T^VF*) is estimated from the 
measurement of the WWH coupling [64]. The branching fractions for the decays to SM 
particles, in particular, bh and 77 would also be lower by an order of magnitude, making 
Higgs searches at the LHC more difficult. 

Therefore, it seems that this model of mixed sneutrino DM has signals at colliders, and 
is testable through direct and indirect detection experiments over much of the parameter 
space. 

5. Conclusions 

The model detailed herein constitutes a viable model of both DM and neutrino mass gen- 
eration in which the DM is comprised of mixed Ihd-rhd sneutrinos in the same flavour mix- 
tures as the neutrino mass eigenstates. We have shown in detail that it is possible that the 
observed neutrino masses and mixings arise as a consequence of supersymmetry breaking 
effects in the sneutrino sector, consistent with all experimental constraints. The prospects 
for indirect detection of the associate sneutrino DM from solar capture are good, owing 
to the large sneutrino-nucleon cross-section afforded by the inelastic splitting of Z-coupled 
states as shown in Figures 5, 6, and 7. For some regions of the allowed parameter space a 
potential 'smoking gun' signature of this model would be monochromatic, non-oscillating, 
neutrinos comprised of just one mass eigenstate which originate from sneutrino annihilation 
in the Sun, see Figure 11. The current generation of direct detection experiments place 
strong limits at low sneutrino masses and at high mass if the inelastic splitting is small, 
and next generation experiments can explore much of the parameter space, see Figures 5, 

6, and 7. There is a prospect of seeing both elastic and inelastic scattering in these exper- 
iments, allowing the size of the sneutrino splitting 5a,L to be directly measured. Despite 
being one of the earliest examples of inelastic DM, limits from solar capture completely 
exclude this model as an explanation of DAMA/LIBRA. The mixed Ihd-rhd sneutrinos can 
have signatures at the LHC, and can change aspects of Higgs physics. Finally, we have 
shown that a simple extension of the MSSM leads to a much richer structure in the DM 
sector linking different areas of beyond-the-Standard-Model physics, namely dark matter 
and neutrino physics. 
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A. Sneutrino Masses and Mixings 



After electroweak symmetry breaking the sneutrino mass matrix has the form: 



Mt 



Mil 



Av sin /3I3 
y XijMjyv sin (3 



Av sin PI3 XijMNV sin P\ 

M|l3 XijMNV sin (3 XijM^ 
XijMNV sin (3 Mfls Afsin/3l3 
^i;sin/3l3 ^'^'^ 



M^l3 



(A.l) 



where the basis is {iy* ,n, iy,h*) and individual terms are defined in Section 2.2. If Ux is 
the matrix that diagonahses the A matrix as; Ux ■ X ■ = Diag{Xi, X2, X^), and Pg is a 
permutation matrix given by; 

/I 0\ 
10 
10 
10 
10 

Vo 1/ 

which rotates 6x6 matrices made up of four diagonal 3x3 blocks down to three 2x2 
matrices along the diagonal, then the product; 



(A.2) 



U 




le -le 



l2^Ux 

l2<S)Ux 



(A.3) 



rotates the sneutrino mass matrix down to the block-diagonal form of six 2x2 matrices 
described in Section 2.2. 



B. Neutrino Masses and Mixings 

From the form of the neutrino mass formula one can see that the neutrino mass matrix is 
diagonalized by the matrix Ux to all orders. Considering just the mixing part of eq.(2.8) 
we have the term: 

12 

Mu,, oc J2 ulaL{r^^ rn^JUc.,j+e (B.l) 
0=1 

L{ra, T^xx) is a diagonal matrix and, as Pq is a permutation matrix, Pg • Diag{6 x 6) • Pg is 
also diagonal. Therefore the only non-diagonal contribution to the neutrino mass matrix 
is from the Ux part of the sneutrino mixing matrix and Ux ■ My ■ u\ must be diagonal. 

B.l Perturbative Mass Formula 

Here we give the full form of eq.(2.10). Defining Xx = Mp^/m,^^ and Ux = M^j^/m^^, then: 




(B.2) 



• (sin^ 0(cos^ — sin26)Li{xxiyx) + cos^ ^(sin^ + sin 29)L2{xx, Ux)) 
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where 

1 /x^lnx^ y'^\ny'^\ 1 / Iny^ 

-^1(2;, y) = —2 2 1 2 ~ 1 2" ~ 1 2 1 + 1 2 vB.3) 

x'^ — \i — x'^ I — y^ J I — y^ \ 1 — 



and 



References 

[1] D. B. Kaplan, A Single explanation for both the baryon and dark matter densities, Phys. Rev. 
Lett. 68 (1992) 741-743. 

[2] D. Hooper, J. March-Russell, and S. M. West, Asymmetric sneutrino dark matter and the 
^b/^DM puzzle, Phys. Lett. B605 (2005) 228-236, [hep-ph/0410114]. 

[3] G. R. Farrar and G. Zaharijas, Dark matter and the baryon asymmetry of the universe, 
hep-ph/0406281. 

[4] R. Kitano and I. Low, Dark matter from baryon asymmetry, Phys. Rev. D71 (2005) 023510, 
[hep-ph/0411133]. 

[5] N. Cosme, L. Lopez Honorez, and M. H. G. Tytgat, Leptogenesis and dark matter related?, 
Phys. Rev. D72 (2005) 043505, [hep-ph/0506320]. 

[6] D. Suematsu, Dark world and baryon asymmetry from a common source, JCAP 0601 (2006) 
026, [astro-ph/0511742]. 

[7] S. Abel and V. Page, Affleck-Dine (pseudo)-Dirac neutrinogenesis, JHEP 05 (2006) 024, 
[hep-ph/0601149]. 

[8] S. Gopalakrishna, A. de Gouvea, and W. Porod, Right-handed sneutrinos as nonthermal dark 
matter, JCAP 0605 (2006) 005, [hep-ph/0602027]. 

[9] J. McDonald, Right-Handed Sneutrino Condensate Cold Dark Matter and the 
Baryon-to- Dark Matter Ratio, JCAP 0701 (2007) 001, [hep-pli/0609126]. 

[10] V. Page, Non-thermal right-handed sneutrino dark matter and the floAi/^b problem, JHEP 
04 (2007) 021, [hep-ph/0701266]. 

[11] F. Deppisch and A. Pilaftsis, Thermal Right-Handed Sneutrino Dark Matter in the Fjj-Term 
Model of Hybrid Inflation, JHEP 10 (2008) 080, [arXiv: 0808. 0490]. 

[12] D. E. Kaplan, M. A. Luty, and K. M. Zurek, Asymmetric Dark Matter, Phys. Rev. D79 
(2009) 115016, [arXiv:0901.4117]. 

[13] G. D. Kribs, T. S. Roy, J. Terning, and K. M. Zurek, Quirky Composite Dark Matter, 
arXiv: 0909. 2034. 

[14] T. Cohen and K. M. Zurek, Leptophilic Dark Matter from the Lepton Asymmetry, 
arXiv: 0909. 2035. 

[15] N. Arkani-Hamed, L. J. Hall, H. Murayama, D. Tucker-Smith, and N. Weiner, Small neutrino 
masses from supersymmetry breaking, Phys. Rev. D64 (2001) 115011, [hep-ph/0006312]. 

[16] N. Arkani-Hamed, L. J. Hall, H. Murayama, D. Tucker-Smith, and N. Weiner, Neutrino 
masses at v**(3/2), hep-ph/0007001. 



- 22 - 



[17] F. Borzumati and Y. Nomura, Low-scale see-saw mechanisms for light neutrinos, Phys. Rev. 
D64 (2001) 053005, [hep-ph/0007018]. 

[18] F. Borzumati, K. Hamaguchi, Y. Nomura, and T. Yanagida, Variations on supersymmetry 
breaking and neutrino spectra, hep-ph/0012118. 

[19] D. Tucker-Smith and N. Weiner, Inelastic dark matter, Phys. Rev. D64 (2001) 043502, 
[hep-ph/0101138]. 

[20] D. Tucker-Smith and N. Weiner, The status of inelastic dark matter, Phys. Rev. D72 (2005) 
063509, [hep-ph/0402065]. 

[21] J. March-Russell and S. M. West, A simple model of neutrino masses from supersymmetry 
breaking, Phys. Lett. B593 (2004) 181-188, [hep-ph/0403067]. 

[22] T. Hambye, J. March-Russell, and S. M. West, TeV scale resonant leptogenesis from 
supersymmetry breaking, JHEP 07 (2004) 070, [hep-ph/0403183]. 

[23] S. M. West, Naturally degenerate right handed neutrinos, Phys. Rev. D71 (2005) 013004, 
[hep-ph/0408318]. 

[24] DAMA Collaboration, R. Bernabei et al. First results from DAMA/LIBRA and the 

combined results with DAMA/Nal, Eur. Phys. J. C56 (2008) 333-355, [arXiv: 0804. 2741]. 

[25] Y. Cui, D. E. Morrissey, D. Poland, and L. Randall, Candidates for Inelastic Dark Matter, 
JHEP 05 (2009) 076, [arXiv: 0901 . 0557]. 

[26] S. Nussinov, L.-T. Wang, and I. Yavin, Capture of Inelastic Dark Matter in the Sun, 
arXiv: 0905. 1333. 

[27] A. Menon, R. Morris, A. Pierce, and N. Weiner, Capture and Indirect Detection of Inelastic 
Dark Matter, arXiv : 0905 . 1847. 

[28] L. J. Hall, T. Moroi, and H. Murayama, Sneutrino cold dark matter with lepton-number 
violation, Phys. Lett. B424 (1998) 305-312, [hep-ph/9712515]. 

[29] S. Kolb, M. Hirsch, H. V. Klapdor-Kleingrothaus, and O. Panella, Collider signatures of 
sneutrino cold dark matter, Phys. Lett. B478 (2000) 262-268, [hep-ph/9910542]. 

[30] T. Asaka, K. Ishiwata, and T. Moroi, Right-handed sneutrino as cold dark matter, Phys. Rev. 
D73 (2006) 051301, [hep-ph/0512118]. 

[31] T. Asaka, K. Ishiwata, and T. Moroi, Right-handed sneutrino as cold dark matter of the 
universe, Phys. Rev. D75 (2007) 065001, [hep-ph/0612211]. 

[32] H.-S. Lee, K. T. Matchev, and S. Nasri, Revival of the thermal sneutrino dark matter, Phys. 
Rev. D76 (2007) 041302, [hep-ph/0702223]. 

[33] C. Arina and N. Fornengo, Sneutrino cold dark matter, a new analysis: Relic abundance and 
detection rates, JHEP 11 (2007) 029, [arXiv: 0709. 4477]. 

[34] Z. Thomas, D. Tucker-Smith, and N. Weiner, Mixed Sneutrinos, Dark Matter and the LHC, 
Phys. Rev. D77 (2008) 115015, [arXiv : 0712 . 4146]. 

[35] C. Arina, Sneutrino cold dark matter in extended MSSM models, arXiv : 0805 . 1991. 

[36] C. Arina, F. Bazzocchi, N. Fornengo, J. C. Romao, and J. W. F. Valle, Minimal supergravity 
sneutrino dark matter and inverse seesaw neutrino masses, Phys. Rev. Lett. 101 (2008) 
161802, [arXiv: 0806. 3225]. 



- 23 - 



[37] D. G. Cerdeno, C. Munoz, and O. Seto, Right-handed sneutrino as thermal dark matter, 
Phys. Rev. D79 (2009) 023510, [arXiv : 0807 . 3029]. 

[38] R. AUahverdi, B. Dutta, K. Richardson-McDaniel, and Y. Santoso, Sneutrino Dark Matter 
and the Observed Anomalies in Cosmic Rays, Phys. Lett. B677 (2009) 172-178, 
[arXiv: 0902. 3463]. 

[39] D. G. Cerdeno and O. Seto, Right-handed sneutrino dark matter in the NMSSM, JCAP 0908 
(2009) 032, [arXiv: 0903. 4677]. 

[40] D. A. Demir, L. L. Everett, M. Frank, L. Selbuz, and I. Turan, Sneutrino Dark Matter: 
Symmetry Protection and Cosmic Ray Anomalies, arXiv : 0906 . 3540. 

[41] R. AUahverdi, S. Bornhauser, B. Dutta, and K. Richardson-McDaniel, Prospects for Indirect 
Detection of Sneutrino Dark Matter with IceCube, Phys. Rev. D80 (2009) 055026, 
[arXiv: 0907. I486]. 

[42] D. G. Cerdeno, Thermal right-handed sneutrino dark matter with a singlet Higgs, AIP Conf. 
Proc. 1115 (2009) 163-168. 

[43] R. AUahverdi, Sneutrino dark matter in light of PAMELA, arXiv : 0909 . 5643. 

[44] A. Kumar, D. Tucker-Smith, and N. Weiner, Neutrino Mass, Sneutrino Dark Matter and 
Signals of Lepton Flavor Violation in the MRSSM, arXiv: 0910. 2475. 

[45] WMAP CoUaboration, E. Komatsu et ai, Five-Year Wilkinson Microwave Anisotropy 
Probe (WMAP) ObservationstCosmological Interpretation, Astrophys. J. Suppl. 180 (2009) 
330-376, [arXiv: 0803. 0547]. 

[46] W. Hu and J. Silk, Thermalization constraints and spectral distortions for massive unstable 
relic particles, Phys. Rev. Lett. 70 (1993) 2661-2664. 

[47] L. J. HaU, K. Jedamzik, J. March-RusseU, and S. M. West, Freeze-In Production of FIMP 
Dark Matter, arXiv:0911 . 1120. 

[48] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, Dark matter direct detection rate in 
a generic model with micrOMECAs2.1, Comput. Phys. Commun. 180 (2009) 747-767, 
[arXiv: 0803. 2360]. 

[49] A. Semenov, LanHEP - a package for the automatic generation of Feynman rules in field 
theory. Version 3.0, Comput. Phys. Commun. 180 (2009) 431-454, [arXiv : 0805 . 0555]. 

[50] Super-Kamiokande Collaboration, S. Desai et ai. Search for dark matter WIMPs using 
upward through-going muons in Super-Kamiokande, Phys. Rev. D70 (2004) 083523, 
[hep-ex/0404025]. 

[51] J. March-RusseU, C. McCabe, and M. McCuUough, Inelastic Dark Matter, Non-Standard 
Halos and the DAMA/LIBRA Results, JHEP 05 (2009) 071, [arXiv: 0812 . 1931]. 

[52] CDMS Collaboration, Z. Ahmed et ai. Search for Weakly Interacting Massive Particles with 
the First Five-Tower Data from the Cryogenic Dark Matter Search at the Soudan 
Underground Laboratory, Phys. Rev. Lett. 102 (2009) 011301, [arXiv: 0802. 3530]. 

[53] XENON CoUaboration, J. Angle et ai. First Results from the XENON 10 Dark Matter 
Experiment at the Cran Sasso National Laboratory, Phys. Rev. Lett. 100 (2008) 021303, 
[arXiv: 0706. 0039]. 



- 24 - 



[54] S. Chang, G. D. Kribs, D. Tucker- Smith, and N. Weiner, Inelastic Dark Matter in Light of 
DAMA/LIBRA, Phys. Rev. D79 (2009) 043513, [arXiv : 0807 . 2250]. 

[55] D. B. Chnc, W. Ooi, and H. Wang, Constraints on Inelastic Dark Matter Signal using 
ZEPLIN- II Results, arXiv:0906.4119. 

[56] K. Schmidt-Hoberg and M. W. Winkler, Improved Constraints on Inelastic Dark Matter, 
JCAP 0909 (2009) 010, [arXiv: 0907. 3940]. 

[57] A. Strumia and F. Vissani, Neutrino masses and mixings and., tLep-ph/0606054. 

[58] H. Yuksel, S. Horiuchi, J. F. Beacom, and S. Ando, Neutrino Constraints on the Dark Matter 
Total Annihilation Cross Section, Phys. Rev. D76 (2007) 123506, [arXiv : 0707 . 0196]. 

[59] F. Halzen and D. Hooper, The Indirect Search for Dark Matter with IceCube, 
arXiv: 0910. 4513. 

[60] M. CireUi et ai. Spectra of neutrinos from dark matter annihilations, Nucl. Phys. B727 
(2005) 99-138, [hep-ph/0506298]. 

[61] H. Kraus et ai, EURECA: The European future of dark matter searches with cryogenic 
detectors, Nucl. Phys. Proc. Suppl. 173 (2007) 168-171. 

[62] E. Aprile, L. Baudis, and f. t. X. Collaboration, Status and Sensitivity Projections for the 
XENONIOO Dark Matter Experiment, arXiv: 0902. 4253. 

[63] L. W. Kastens, S. B. Cahn, A. Manzur, and D. N. McKinsey, Calibration of a liquid xenon 
detector with kr-83m. Physical Review C 80 (2009) 045809. 

[64] H. Yamamoto, Physics at International Linear Collider (LLC), J. Phys. Soc. Jap. 76 (2007) 
111014, [arXiv: 0709. 0899]. 



- 25 - 



